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Abstract 

We report major advances in the research program initiated in "Moment-Based Evidence for 
Simple Rational- Valued Hilbert-Schmidt Generic 2x2 Separability Probabilities" (J. Phys. A, 
45, 095305 [2012]). A function P(a), incorporating a family of six hypergeometric functions, all 
with argument g| = (|) 3 , is obtained. It reproduces a series, a = |,2, ... ,32, of sixty-four 
conjectured Hilbert-Schmidt rational- valued generic 2x2 separability probabilities. These exact ra- 
tios are put forth on the basis of systematic, high-accuracy probability-distribution-reconstruction 
computations, employing 7,501 determinantal moments of partially transposed 4x4 density matri- 
ces. A lengthy expression for P(a) containing six generalized hypergeometric functions is initially 
obtained-making use of the FindSequenceFunction command of Mathematica. A remarkably suc- 
cinct reexpression for P(a) is then found, by Qing-Hu Hou and colleagues, using Zeilberger's 
algorithm ("creative telescoping"), For generic (9-dimensional) two-rebit systems, P{\) = §§, (15- 
dimensional) two-qubit, -P(l) = (a value that had been proposed in J. Phys. A, 40, 14279 
[2007] and supported in Intl. J. Mod. Phys. P>, 26, 1250054 [2012]) and (27-dimensional) two- 
quater(nionic)bit systems, P{2) = J^. 
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I. INTRODUCTION 

The predecessor paper pQ -addressing the relatively long-standing 2x2 separability 
probability question [2"HTT] (cf. [T2T[T4] )-consisted largely of two sets of analyses. (In 
the foundational study [2J, "three main reasons of importance" -philosophical, practical 
and physical-were given for studying this question.) The first set of analyses in pQ was 
concerned with establishing formulas for the bivariate determinantal product moments 
(jp PT \ n \p\ k ^ j k, n = 0,1,2,3,..., with respect to Hilbert-Schmidt (Euclidean/flat) mea- 
sure JT5] [ini sec. 14.3], of generic (9-dimensional) two-rebit and (15-dimensional) two-qubit 
density matrices (p). Here p PT denotes the partial transpose of the 4x4 density matrix 
p. Nonnegativity of the determinant |p PT | is both a necessary and sufficient condition for 
separability in this 2x2 setting [T7] . 



In the second set of analyses in [T], the univariate determinantal moments (\p 



polynomial-based probability distribution reconstruction algorithm of Provost [TBI sec. 2] 
(cf. [19J). This yielded estimates of the desired separability probabilities. (The reconstructed 
probability distributions based on |p PT | are defined over the interval |p PT | G [— ^, ^g], while 
the associated separability probabilities are the cumulative probabilities of these distribu- 
tions over the nonnegative subinterval \p PT \ G [0, ^g]. We note that for the fully mixed 
(classical) state, |p PT | = ^r, while for a maximally entangled state, such as a Bell state, 

IP PT ! = 

A highly-intriguing aspect of the (not yet rigorously established) determinantal moment 
formulas obtained (by C. Dunkl) in (TJ App.D.4] was that both the two-rebit and two-qubit 
cases could be encompassed by a single formula, with a Dyson-index-like parameter a [20] 
serving to distinguish the two cases. The value a = \ corresponded to the two-rebit case 
and a = 1 to the two-qubit case. (Let us note that the results of the formula for a = 2 
and n = 1 and 2 have recently been confirmed computationally by Dunkl using the " Moore 
determinant" (quasideterminant) [21, [22] of 4 x 4 quaternionic density matrices. However, 
tentative efforts of ours to verify the a = 4 [conjecturally, octonionic [23J, problematical] 
case, have not proved successful.) 

When the probability-distribution-reconstruction algorithm [18] was applied in p] to the 
two-rebit case (a = |), employing the first 3,310 moments of |p PT |, a (lower-bound) estimate 





formulas, served as input to a Legendre- 
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that was 0.999955 times as large as §f » 0.453120 was obtained (cf. |24[ p. 6]). 

Analogously, in the two-qubit case (a = 1), using 2,415 moments, an estimate that 
was 0.999997066 times as large as J| ~ 0.242424 was derived. This constitutes an ap- 
pealingly simple rational value that had previously been conjectured in a quite different 
(non-moment-based) form of analysis, in which "separability functions" had been the main 
tool employed [11]. (Note, however, that the the two-rebit separability probability conjec- 
ture of j^, somewhat secondarily advanced in |llj . has now been discarded in favor of ||.) 
Let us, importantly, note that in an extensive Monte Carlo analysis, Zhou, Chern, Fei and 
Joynt obtained an estimate for this two-bit separability probability of 0.2424 ± 0.0002 [231 
eq. (B7)], thus, fully consistent with our exact proposal. 

Further, the determinantal moment formulas advanced in [1] were then applied with a 
set equal to 2. This appears-as the indicated recent (Moore determinant) computations of 
Dunkl show-to correspond to the generic 2 7- dimensional set of quaternionic density matrices 
[2SH2Z]- Quite remarkably, a separability probability estimate, based on 2,325 moments, that 
was 0.999999987 times as large as J|| ~ 0.0804954 was found. (In line with this set of three 
results, the paper p] was entitled, "Moment-Based Evidence for Simple Rational- Valued 
Hilbert-Schmidt Generic 2x2 Separability Probabilities".) 

In the present study, we extend these three (individually-conducted) moment-based anal- 
yses in a more systematic, thorough manner, jointly embracing the sixty-four integral and 
half- integral values a = |, 1, §, 2, . . . , 32. We do this by accelerating, for our specific pur- 
poses, the Mathematica probability-distribution-reconstruction program of Provost [IB] , in 
a number of ways. Most significantly, we make use of the three-term recurrence relations for 
the Legendre polynomials. Doing so obviates the need to compute each successive higher- 
degree Legendre polynomial ab initio. 

In this manner, we were able to obtain-using exact computer arithmetic throughout- 
" generalized" separability probability estimates based on 7,501 moments for a = 
|, 1, |, . . . , 32. In Fig. [l]we plot the logarithms of the resultant sixty-four separability prob- 
ability estimates (cf. [H Fig. 8]), which fall close to the line -0.9464181889a. In Fig. [2] we 
show the residuals from this linear fit. 

In Fig. [3] we present a hypergeometric-f unction-based formula, together with striking sup- 
porting evidence for it, that appears to succeed in uncovering the functional relation (P(a)) 
underlying the entirely of these sixty-four generalized separability probabilities. Further, in 
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FIG. 1: Logarithms of generalized separability probability estimates, based on 7,501 Hilbert- 
Schmidt moments of Ip^l, as a function of the Dyson-index- like parameter a 
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FIG. 2: Residuals from linear fit to logarithms of generalized separability probability estimates 

(j3J), and the immediately preceding text, we list a number of remarkable values yielded by 
this hypergeometric formula for values of a other than the basic sixty-four (half-integral and 
integral) values from which we have started. 

Then, we are able to report-with the assistance of Qing-Hu Hou and colleagues-of a 
remarkable condensation of the lengthy expression presented in Fig. [3j 
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FIG. 3: Hypergeometric formula P(a) for Hilbert-Schmidt generic 2x2 generalized separability- 
probabilities and evidence that it reproduces the basic three (real [a = 5], complex [a = 1] and 
quaternionic [a = 2]) conjectures of ||, J| and 



II. RESULTS 



A. The three basic conjectures revisited 

1. a = 2~the two-rebit case 

In [1] , a lower-bound estimate of the two-rebit separability probability was obtained, with 
the use of the first 3,310 moments of |p PT |- It was 0.999955 times as large as §§ » 0.453120. 
With the use, now, of 7,501 moments, the figure increases to 0.999989567. This outcome, 
thus, fortifies our previous conjecture. 

2. a = 1-the two-qubit case 

In [T], a lower-bound estimate of the two-qubit separability probability was obtained, 
with the use of the first 2,415 moments of |p PT |, that was 0.999997066 times as large as 

w 0.242424 (cf. [251 eq. (B7)]). Employing 7,501 moments, this figure increases to 
0.99999986. 

3. a = 2-the quaternionic case 

In [1] , a lower-bound estimate of the (presumptive) quaternionic separability probability 
was obtained that was 0.999999987 times as large as g « 0.0804954, using the first 2,325 
moments of |p PT |. Based on 7,501 moments, this figure increases, quite remarkably still, to 
0.999999999936. 

B. Generalized separability probability formula 

A principal motivation in undertaking the analyses reported here-in addition, to further 
scrutinizing the three specific conjectures reported in [T]-was to uncover the functional 
relation underlying the curve in Fig. [T] (and/or its original non- logarithmic counterpart). 

Preliminarily, let us note that the zeroth-order approximation (being independent of the 
particular value of a) provided by the Provost probability-distribution-reconstruction al- 
gorithm is simply the uniform distribution over the interval [—-h, ocr]- ^ ne corresponding 
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zeroth-order separability probability estimate is the cumulative probability of this distribu- 
tion over the nonnegative subinterval [0, ^g], that is, ^-/(jq + ^b) = jj ~ 0.0588235. So, 



it certainly appears that speedier convergence (sec. II A ) of the algorithm occurs for separa- 
bility probabilities, the true values of which are initially close to ^ (such as ^ ~ 0.0804954 
in the quaternionic case). Convergence also markedly increases as a increases. 

It appeared, numerically, that the generalized separability probabilities for integral and 
half- integral values of a were rational values (not only ||, for the three specific values 

a = h, 1, 2 of original focus). With various computational tools and search strategies based 
upon emerging mathematical properties, we were able to advance additional, seemingly 
plausible conjectures as to the exact values for a = 3,4, ...,32, as well. (We inserted 
many of our high-precision numerical estimates into the search box on the Wolfram Alpha 
website-which then indicated likely candidates for corresponding rational values.) 

We fed this sequence of thirty-two conjectured rational numbers into the FindSequence- 
Function command of Mathematica. (This command "attempts to find a simple function 
that yields the sequence a, when given successive integer arguments," but apparently can 
succeed with rational arguments, as well.) To our considerable satisfaction, this produced 
a generating formula (incorporating a diversity of hypergeometric functions of the P -F p _i 
type, p = 7, ... ,11, all with argument z = || — (|) 3 ) f° r the sequence. (Let us note that 
z -\ = is the "residual entropy for square ice" [281 p. 412] (cf. [291 eqs.[(27), (28)]. 
An analogous appearance of || occurs in a hypergeometric ["Ramanujan-like"] summation 



for 1^- [30J of Guillera. In a related email he wrote: "Some values of "z", for example z 
= 27 / 64 appear frequently in hypergeometric identities. I think it is due to a "modular" 
or a "modular-like" origin."). In fact, the Mathematica command succeeds using only the 
first twenty-eight conjectured rational numbers, but no fewer-so it seems fortunate, our 
computations were so extensive.) 

However, the formula produced by the Mathematica command was quite cumbersome 
in nature (extending over several pages of output). With its use, nevertheless, we were 
able to convincingly generate rational values for /iaZ/-integral a (including the two-rebit 
|| conjecture), also fitting our corresponding half-integral thirty- two numerical estimates 
exceedingly well. (Let us strongly emphasize that the hypergeometric-based formula was 
generated using only the integral values of a. The process was fully reversible, and we could 
first employ the half-integral results to generate the formula-which then-seemingly perfectly 
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fitted the integral values.) 

At this point, for illustrative purposes, let us list the first ten half- integral and ten in- 
tegral rational values (generalized separability probabilities), along with their approximate 
numerical values. 



a = | || 0.453125 a = 1 ^ 0.242424 

a = | ^jga 0.137562 a = 2 ^ 0.0804954 

« = I iHil 0.0480083 a=3jS 0.0290081 

« = I 2i990^5555 2 0-0 1 76 9 48 a = 4 ^ 0.0108722 

« = 1 90on99So 9 92 0.00671852 a = 5 0.00417101 

_ 11 71925602948804923 D 009 ^QQ4 rv = fi 179808469 OOOIf^lQ 

" 2 27670116110564327424 "-uu^^** « u 110638410169 u - ,J,J - LUZ " JJ - t ' 

13 3387374833367307236269 n DDI 01 RQ rv = 7 191151001 n OOOfUO'WQ 

u 2 3324546003940230230441984 u - uuiuios « ' 298529164591 U - UUUU ^ U ' JUJ 

15 124792688228667229196729 n 00040^997 rv = 8 1331199762 r, 0009^01 

2 309485009821345068724781056 u.uuu^uo^z, 1 « ° 5232880523393 u.uuuz^ojx 

17 407557367133399293946182513 0001 fi07^ rv = Q 74195568677 0001 01 79Q 

2 2535301200456458802993406410752 u.uuuxuu i ul j 729345064647247 ,J - ,J,J,J - L,J - L 1 



19 1338799759394288468677657208071 n 0000R44R0Q rv = 1 730710456538 n 00004080^0 
u 2 20769187434139310514121985316880384 u.uuuuu<±<±uu» « lu 17868447453498669 u - UUUU4:UO » d » 

(1) 

To simplify the cumbersome (several-page) output yielded by the Mathematica FindSe- 
quenceFunction command, we employed certain of the " contiguous rules" for hypergeometric 
functions listed by C. Krattenthaler in his package HYP [3T] (cf. [32]). Multiple applica- 
tions of the rules C14 and C18 there, together with certain gamma function simplifications 
suggested by C. Dunkl, led to the rather more compact formula displayed in Fig. |3j This 
formula incorporates a six-member family {k = 1, . . . , 6) of ^Fq hypergeometric functions, 
differing only in the first upper index k, 

„ /, 2 3 4 5 7 6 13 3 17 19 21 27 

7 F 6 ^k,a+-,a + -,a+-,a + -,a + -,a+pa + —,a + -,a+—,a+—,a + 2,a+— ] — 

(2) 

(The reader will note interesting sequences of upper and lower parameters (cf. [33]).) We 
are only able to, in general, evaluate the formula numerically, but then to arbitrarily high 
(hundreds, if not thousand-digit) precision, giving us strong confidence in the validity of the 
exact generalized separability probabilities that we advance. 

Let us now apply the formula (Fig. [3]) to values of a other than the basic sixty-four. For 
a = 0, the formula yields-as would be expected-the "classical separability probability" of 1. 
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Further, proceeding in a purely formal manner (since there appears to be no corresponding 
genuine probability distribution over [— i, ^)), for the negative value a = — |, the formula 
yields |. For a = — |, it gives -2. Remarkably still, for a = |, the result is clearly (to 
one thousand decimal places) equal to 2 — 



34 



, _ - .. ., ^ 0.6486993992, where 

21agm(l,V2j 21V27T 3 / 2 ' 

the arithmetic-geometric mean of 1 and \f2 is indicated. (The reciprocal of this mean is 

9689rf — i 

Gauss's constant.) For a — § , the result equals 2 ^ ,{\ ~ 0.3279684732, while for 



128 



+ 2 = 2 + 



21tt3/ 2 



4420^(1) 

7.087249321. For a 



a = — f , we have ; — — 

4' 21agm(l,y2) 
288927rf — 1 

^ 344080,1-2 ~ 0.36424897456. Results are presented in the table 



, the outcome is 



a 

_3 

4 

_ 2 
3 

_1 

2 

_ 1 

3 

_ 1 

4 

1 
1 

1 

3 

2 
3 



P(a) 

32V2r(|^ 

217T3/2' 
8tt 

" Vsr(i) : 

2 

2+ 3r(|) 3 



2 + 
2 - 



2 - 
2 



3 



2 
2 



4tt 2 



21 V / 2W 2 

_ 459^^ 

288927r(i)' 
3440807T 2 
9689r(|) 

4420v^r(|^ 



value \ 

7.08725 
1.24527 

0.666667 

3.461 
2 

0.648699 
0.572443 

0.364249 
0.327968 



(3) 



(Let us note that the term 



3r 



4tt 2 



1.46099848 present in the result for a = — \ is "Baxter's 



four-coloring constant" for a triangular lattice [281 P- 413].) Also, for a 



1 :,. " i 

-1, we have f. 
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For a — — §, the result is |. 



III. CONCISE FORMULATION OF HYPERGEOMETRIC EXPRESSION 



The concise formulation of P(ot) takes the form (Fig. |4 



P(a) = Z°° =0 f(a + i) 



(4) 



where 



f{a) = P{a) - P{a + 1 = VTuvo Z quv* ^ i3x 

3r(a + l)T(2a + 3)T(5a + — i 



l3^ ' 

2 



(5) 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 



a 



FIG. 4: Generalized two-qubit separability probability function P(a), with P(0) = l,P(h) = 
29 ,P(1) = Jr,P(2) = Jpj for generic classical four-level (a = 0), two-rebit (a = i), two-qubit 



64' 



(a = 1) and two-quaterbit (a = 2) systems, respectively. 



and 



q(a) = 185000a 5 + 779750a 4 + 1289125a 3 + 1042015a 2 + 410694a + 63000 = (6) 
a(5a(25a(2a(740a + 3119) + 10313) + 208403) + 410694) + 63000. 



We had previously ourselves been unable to find an equivalent form of P{a) with fewer 
than six hypergeometric functions [3U Fig. 3] . Qing-Hu Hou and colleagues of the Center for 
Combinatorics of Nankai University, however, were able to obtain the remarkably succinct 
and clearly correct results Q- ^ -which they communicated to us in a few brief, informal 
e-mail messages. (They provided two Maple worksheets they employed [Figs. [5] and |6j.) 
They, first, observed that the hypergeometric-based formula for P(a) could be expressed as 
an infinite summation. Letting Pi(a) be the /-th such summand, application of Zeilberger's 
algorithm [35] (a method for producing combinatorial identities, also known as "creative 
telescoping") yielded that 



P t (a) - P t {a + 1) = -Pi +1 {a) + P,(a). 



(7) 



(The package APCI-available at http://www.combinatorics.net.cn/homepage/hou/-was 
employed. In a different quantum-information context, Datta employed the algorithm to 
ascertain that no closed form exists for a certain series, "retarding" the evaluation of the 
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ratio "of the negativity of random pure states to the maximal negativity for Haar-distributed 
states of n qubits [SBl App. A].) Summing over I from to oo, Hou and colleagues found 
that 

P(a) - P(a + 1) = P (a). (8) 

Letting f(a) = Po(ot), the concise summation formula ^ is obtained. (C. Krattenthaler 
indicated-and Hou agreed-that these results might equally well be derived without recourse 
to Zeilberger's algorithm. Also, a referee expressed puzzlement at the peculiar form of 
eq. ([7]). This appears to be an artifact arising from the particular manner in which the 
algorithm is applied in the proving of hypergeometric identities.) 
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> with(APCI) ; 

[AbelZ, Disset, ExtZeil, Gosper, Zeil, hyper _simp, hyper term, poch, qExtJZeil, qGosper, (1) 
qZeil, qhyper _simp, qhyperterm, qpoch] 

> t:=hyperterm( [k,a+2/5,a+3/5,a+4/5,a+5/6,a+7/6,a+6/5] , [a+13/10, 
a+3/2,a+17/10,a+19/10,a+2,a+21/10] ,27/64,1) ; 

pochhammer (k, I) pochhammer i^a + / j pochhammer fa + pochhammer fa (2) 

~*~ T' '] P oc ' 1 ' iarnmer [ a ~*~ '] pochhammer fa + / I pochhammer fa + 



? := 



/ 



27_ 
64 



pochhammer f a + -j^-, / j pochhammer f a + / 



17 ^ f 19 "\ f 21 

+ -j^-, / pochhammer a + —, I pochhammer (a + 2, 1) pochhammer a + —, 

,)„) 

> PP:=(a*(5*a*(25*a*(2*a*(740*a-581)+161)+628)+39)-54)*subs(k=l,t)+ 
(5*a*(25*a*(8*a*(925*a-2431) +22255) -312019) +347274 )*subs(k=2,t) 
+10* ( ( 25*a*(4*a*(3700*a-12843) +66227) -769797 )*subs(k=3,t) +75* ( 
(8*a*(1850*a-6131)+44133)*subs(k=4,t) + 8*( (3700*a-7981) *subs(k= 
5,t) + 3700*subs(k=6,t) ) ) ): 

> Pl:=4"(-2*a-3) *GAMMA(3*a+5/2)*GAMMA(5*a+2) /3/GAMMA(a+l)/GAMMA(2* 
a+3) /GAMMA(5*a+13/2) *PP: 

> re:=Zeil(Pl,a,l, "cert" ) ; 

re:= [S(a) -S(a + 1) =0, -1] (3) 



> f :=- factor (hyper_s imp ( subs(l=0,Pl*re[2] ) )); 

1 

a\ (2a + 2)! f 5 a + 



/:= j - ( (185000 a + 779750 a 4 + 1289125 a + 1042015 a 2 (4) 



2 



410694 a + 63000) 2~ 4 17 6 (3a+|-|!(5a+l)! 



> add(eval£ (hyper_simp(subs (a=l/4+i, f))), i=0..30); 

0.6486993992 (5) 

> add(evalf (hyper_simp(subs(a=-l+i, f))), i=0..30); 

-2.000000001 (6) 

> add(evalf (hyper_s imp (subs (a=i, f))), i=0..30); 

0.9999999999 (7) 
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FIG. 5: First Maple worksheet of Hou used in deriving concise form of hypergeometric formula 



> with(APCI); 

[AbelZ, Disset, Ext_Zeil, Gosper, Zeil, hyper _simp, hyper term, poch, qExtJZeil, qGosper, 
qZeil, qhyper_simp, qhyperterm, qpoch ] 

"> t:=hyperterm( [k,a+2/5,a+3/5,a+4/5,a+5/6,a+7/6,a+6/5] , [a+13/10, 
a+3/2,a+17/10,a+19/10,a+2,a+21/10] ,27/64,1) ; 



t := 



2 ^ ( 3 



pochhammer (k,l) pochhammer + —, l\ pochhammer yi + — , / J pochhammer la 
f / j pochhammer (a + / I pochhammer fa + / I pochhammer fa + 



/ 



13 ,\ ( 3 



[ pochhammer | a + — , / pochhammer \a + — , / | pochhammer | a 



' 27 
,64 

17 ^ 19 ^ (21 

-j^-, / pochhammer a + —, I pochhammer {a + 2,1) pochhammer a + 



I l\ 



> PP:=(a*(5*a*(25*a*(2*a*(740*a-581)+161)+628)+39)-54)*subs(k=l,t) 
(5*a*(25*a*(8*a*(925*a-2431) +22255) -312019) +347274 )*subs(k=2,t) 
+10* ( ( 25*a*(4*a*(3700*a-12843) +66227) -769797 )*subs(k=3,t) +75* ( 
(8*a*(1850*a-6131)+44133)*subs(k=4,t) + 8*( (3700*a-7981)*subs(k 
5,t) + 3700*subs(k=6,t) ) ) ): 

"> Pl:=4"(-2*a-3) *GAMMA(3*a+5/2) *GAMMA(5*a+2)/3/GAMMA(a+l)/GAMMA(2* 
a+3) /GAMMA(5*a+13/2 ) *PP: 



> re:=Zeil(Pl,a,l,"cert" ) ; 

re:= [S(a) -S(a + 1) = 0, -1] 

> zz :=-f actor (hyper_simp( subs(l=0,Pl*re[2] ) )) 



1 1 



3 ( 11 

a\ (2a + 2)! 5 a + — | ! 



(185000 a 5 + 779750 a + 1289125 a 



+ 1042015 a 2 + 410694 a + 63000) 2~ 4a 6 (3 a + j j ! (5 a + 1 ) ! 

The ratio of two consecutive terms 

> rl:=hyper_simp( subs (1=1+1, PI) /PI) ; 

rl := (3 (3769584 + 10406099 a + 10406099 /+ 1704750 a 4 + 11437890 a 2 + 6258125 a 
+ 185000 a + 925000 a 4 / + 6819000 a I + 18774375 a 2 / + 22875780 a / 
+ 1850000 a 3 I 2 + 10228500 a 2 I 2 + 18774375 a/ 2 + 11437890 I 2 + 1850000 a 2 / 3 
+ 6819000 a / 3 + 6258125 / 3 + 925000 a I 4 + 1704750 I 4 + 185000 I 5 ) (5 a + 2 
+ 5/) (5a + 3 + 5/ ) (5a + 4 + 5/ ) (6 a + 5 + 6/) (6a + 7 + 6/) (5a + 6 
+ 5 /))/(8 (10 a + 21 + 10/) (a + 2 + / ) (10 a + 19 + 10/) (10 a + 17 + 10/) (2 a 
+ 3 +2 1) (10 a + 13 + 10/) (63000 + 410694 a + 410694/ + 779750 a 4 
+ 1042015 a 2 + 1289125 a + 185000 a + 925000 a I + 3119000 a I + 3867375 a 2 / 
+ 2084030 a / + 1850000 a I 2 + 4678500 a 2 I 2 + 3867375 a I 2 + 1042015 I 2 
+ 1850000 a / 3 + 3119000 a/ 3 + 1289125 / 3 + 925000 a/ 4 + 779750 I 4 + 185000 / 5 ) 



13 



= ) 

> r2:=hyper_simp( subs (a=a+l+l,zz) /subs (a=a+l,zz) ) ; 

r2 := (3 (3769584 + 10406099 a + 10406099 / + 1704750 a + 11437890 a + 6258125 a (6) 

+ 185000 a + 925000 a 4 / + 6819000 a I + 18774375 a 2 I + 22875780 a I 

+ 1850000 a 3 I 2 + 10228500 a I 2 + 18774375 a I 2 + 11437890 I 2 + 1850000 a / 3 

+ 6819000 a / 3 + 6258125 / 3 + 925000 a I 4 + 1704750 I 4 + 185000 / 5 ) (5 a + 2 
+ 5/) (5 a + 3 + 5/ ) (5a + 4 + 5/ ) (6a + 5 + 6/) (6a + 7 + 6/) (5 a + 6 

+ 5/))/ (8 (10 a + 21 + 10/) (a + 2+l) ( 10 a + 19 + 10/) (10 a + 17+ 10/) (2 a 

+ 3 + 2 /) (10 a + 13 + 10 /) (63000 + 410694 a + 410694 / + 779750 a 4 

+ 1042015 a + 1289125 a + 185000 a + 925000 a I + 3119000 a I + 3867375 a I 

+ 2084030 a / + 1850000 a I 2 + 4678500 a I 2 + 3867375 a I 2 + 1042015 I 2 

+ 1850000 a 2 / 3 + 31 19000 a / 3 + 1289125 / 3 + 925000 a I 4 + 779750 / 4 + 185000 / 5 ) 

= ) 

> normal ( rl-r2 ); 

(7) 



14 

FIG. 6: Second Maple worksheet of Hou used in deriving concise form of hypergeometric formula 



We certainly need to indicate, however, that if we do explicitly perform the infinite sum- 
mation indicated in Q, then we revert to a (" nonconcise" ) form of P(a), again containing 
six hypergeometric functions. Further, it appears that we can only evaluate Q numerically- 
but then easily to hundreds and even thousands of digits of precision-giving us extremely 
high confidence in the specific rational-valued Hilbert-Schmidt separability probabilities ad- 
vanced. 

There remain the important problems of formally verifying the formulas for P(a) (as 
well as the underlying determinantal moment formulas in [Tj, employed in the probability- 
distribution reconstruction process), and achieving a better understanding of what these 
results convey regarding the geometry of quantum states [El EH EE]- Further, questions of 
the asymptotic behavior of the formula (a — > oo) and of possible Bures metric [3| 14"! IB"! IT6l 139] 
counterparts to it, are under investigation |40j. Also, we are presently engaged in attempting 
to determine further properties-in addition to the cumulative (separability) probability over 
[0, 2^g]~of the probability distributions of |p PT | over [—jg, ^g], as a function of the Dyson- 
index-like parameter a. As one finding in this direction, it appears that the ^-intercept 
(at which |p PT | = 0) in the presumed quaternionic case (a = 2) is = ^^j* 11 . (The 
Legendre-polynomial-based probability-distribution reconstruction algorithm of Provost [18] 
yielded an estimate 0.99999999742 times as large as ^||^, when implemented with 10,000 
moments. Based also on 10,000 moments, the two-qubit [a = 1] counterpart was estimated 
as 389.994942076.) 

The foundational paper of Zyczkowski, Horodecki, Sanpera and Lewenstein," Volume of 
the set of separable states" [2], did ask for volumes, not specifically probabilities. At least, 
for the two-rebit, two-qubit and two-quaterbit cases, a = |, 1 and 2, we can readily convert 
the corresponding separability probabilities to the separable volumes gjUf^j = 2^¥^7> 

449513064000 ~~ 2 6 -3 6 -5 3 -7 2 -ll 2 -13 all< ^ 3914156909371803494400000 ~~ 2 14 -3 10 -5 5 -7 3 -ll 2 -13-17 2 -19 2 -23 ' usrn § ^ ne 

Hilbert-Schmidt volume formulas of Andai [2HI Thms. 1-3] (cf. [To], [IB"]). The determination 
of separable volumes-as opposed to probabilities-for other values of a than these three 
appears to be rather problematical, however. 

Let us also note that Theorem 2 of [H], in conjunction with the results here, al- 
lows us to immediately obtain the separability probabilities of the generic minimally- 
degenerate/boundary 8-, 14-, and 26-dimensional two-rebit, two-qubit, and two-quaterbit 
states, as one-half (that is, ^ and ^) the separability probabilities of their generic 
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non-degenerate counterparts. 



IV. APPENDIX EXACT VALUES OF DERIVATIVES OF P(a) 
A. Succeeding deriviatives at a = 

The first derivative of P(a) evaluated at (the classical case) a = is -2, while the second 
derivative is 40 — 20£(2) = 40 — ^rf- ~ 7.10132. (The third derivative was computed as 
-43.7454236566749417600.) 



B. First derivatives at a = 1, 2 . . ., et al 



The first derivative of P(a) at a = -\ is -f and at a = \ is ^(917 - 9841og(2)) w 
0.611831, and -2 at a — 0, as previously mentioned. We have also been able to determine ra- 
tional values of P(ot) for a — 1,2, ... , 97. We list the first seven of these. (The Mathematica 
command FindSequenceFunction, however, did not succeed in this instance in generating an 
underlying function for this sequence of 97 rational numbers-although, of course, one can 
be directly obtained from our explicit forms above of P(a)). 

( rv P'( a ) \ 

« -0.285489 
-0.0845257 
-0.0289266 
w -0.0105241 
-0.00395978 
-0.00152186 
-0.000593494 



a 
1 
2 
3 
4 
5 
6 
7 



457380 

3177826243 
37595998440 

3598754002551529 
124409677632540300 

943222153906869801499 ^ 
89625168823088671652880 ~ 

7745868905935978063871447 
1956135029605259737354520400 

163704960709243940550573265691777 , 
107569184582725029279135417408286275 ' 

124555275071579876642057723808475761407 



(9) 



209867628485254931732709294271962333917400 



/ 
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